Derealization and the semiclassical description of molecular rotation 
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We discuss phase-space derealization for the rigid rotator within a semiclassical context by recourse 
to the Husimi distributions of both the linear and the 3D— anisotropic instances. Our treatment 
is based upon the concomitant Fisher information measures. The pertinent Wehrl entropy is also 
investigated in the linear case. 
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I. INTRODUCTION 

We will be concerned here with the semiclassical descrip- 
tion of the rotational dynamics of molecular systems. A 
pioneer effort in this sense is that of Morales et al. in 
Ref. 0, who studied the connection between such dy- 
namics and coherent states 0,0. 

The coherent states formulation is not unique and sev- 
eral authors have developed alternative descriptions for 
molecular quantum systems. We have chosen two tra- 
ditional formulations for coherent states in order to dis- 
cuss the semiclassical dynamics of molecular rotational 
systems. One of them is devoted to two dimensional case 
Q and the other to the three dimensional one Q. 
In Refs. [H, |ij the authors discuss, among other things, 
the advantages of the Husimi distribution for the inter- 
pretation of the electronic structure in hydrogen and ni- 
trogen atoms. They suggested that their work may be 
extended to the molecular instance. Following this sug- 
gestion we address here the simplest applicable model, 
i.e., the rigid rotator. Its usefulness for describing di- 
atomic molecules is well-known 0. 

Derealization is an "energetically favorable" process, 
since it distributes the wave of function over a vol- 
ume greater than the size of the sample. Thus, the 
net energy of the molecule is lowered, which results in 
resonance-stabilization. The celebrated Fisher informa- 
tion I [1, 0, E3 can De related to the derealization 
measure as follows: if we take a wave package with a 
standard deviation <j the Fisher information is given by 
I > 1/er 2 , thus, in this sense, we realize that the Fisher 
information is a quite sensitive indicator of the dereal- 
ization of the wave package H 13113. 
The rigid rotator is a system of a single particle whose 
quantum spectrum of energy is exactly known. There- 



fore, the study of typical thermodynamic properties can 
be analytically derived^!]]. Applications lead to the 
treatment of important aspects of molecular systems [l2|. 
The paper is organized as follows: In Section [H] we ex- 
plore the linear rigid rotator. We write the probability of 
finding a quantum state in a coherent state that is used 
to obtain an explicit expression for the (i) Husimi distri- 
bution, (ii) Wehrl entropy, and (iii) Fisher Information. 
These results are of help in Section ITTT1 where we discuss 
a model for the three dimensional rigid rotator. Finally, 
in Section ll VI we draw some conclusions. 



II. LINEAR RIGID ROTATOR 

We start the present study by exploring a simple model, 
the linear rigid rotator (LRR), based on the excellent 
discussion concerning the coherent states for angular mo- 
menta given in Ref. 14]. The LRR-hamiltonian writes 



H = 



L 1 



21* 



(1) 



where L 2 = Ll 



T 2 



is the angular momentum operator 



and I x and I y are the associated moments of inertia. We 
have assumed that I xy = I x = I y . Calling \I, K) the set 
of iJ-eigenstates, we recall that they verify the relations 

L 2 \IK) = I{I+l)\IK) 

L Z \IK) = K\IK), (2) 

with I = 0, 1, 2 ... , for —I < K < J, the eigenstates' 



2 



energy spectrum being given by 



si 



i{i + i)h 2 



21 



(3) 



■HJ 



where for simplicity we take h = 1. 

We construct the lineal coherent states for the rigid ro- 
tator using Schwinger's oscillator model of angular mo- 
mentum 0,0] as 



where n + = I + K and n_ = I — K. Therefore, the 
probability of observing the state \IK) in the coherent 
state \ziz 2 ) is of the form 



\{IK\zi*2)t 



2n + 2n_ 
2 z l z 2 



Husimi distribution 



(15) 



\IK) = \ +J V ~ ; 

y/(I + K)\(I-K)\ 



10), 



(4) 



with a+ and a_ the pertinent creation and annihilation 
operators, respectively, and |0) = |0, 0) the vacuum state. 
The states \IK) are orthogonal and satisfy the closure 
relation, i.e., 



(IK\IK) = 6 fj 6 K . K , 



£ £ \IK){IK\ = t. 

I=0K=-I 



(5) 



(G) 



Since we deal with two degrees of freedom the ensuing co- 
herent states are of the tensorial product form (involving 
\zi) and \z 2 )) 



where 



\z\z 2 ) = \zi) ® \z 2 ), 
a+\ziz 2 ) = Zi\ziZ2), 



a-\zxz 2 ) = z 2 \z x z 2 ). 
Therefore, the coherent state \z\Z 2 ) writes [Ti) 



with 



\ziz 2 ) 
\zi) 

l*2> 



1*1 -t -t 

+ e Zia +e Z2a -\0), 



^e zia +|0), 

t 



1 a 



-|0>. 



We have introduced the convenient notation 



\z? = N 2 + N 2 - 



(7) 
(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



Using Eqs. (0J and (|10|l we easily calculate llK\z\z 2 ) 
and, after a bit of algebra, find 



Following the procedure developed by Anderson et 
al. Ilal. we can readily calculate the Husimi distribu- 
tion ml, which is defined as 



H(zi,z 2 ) = (zi,z 2 \p\zi,z 2 ), 
where the density operator is 



p = Z 2 p exp ( 



(16) 



(17) 



and j3 = 1/fceT, ks is the Boltzmann's constant and T 
the temperature. The form of the rotational partition 
function Z 2 u is given in Ref. Q 



Ji) =^(2/+l) e -^+ 1 )f : 

7=0 



(18) 



with 9 = h 2 1 (2I xy kB)- In the present context, speaking 
of the "trace operation" entails performing the sum Tr = 

EDO IT 1 / 

1=0 1^K=-I ■ 

Inserting the closure relation into Eq. (|16|l . and using 
Eq. Q15|). we finally get our Husimi distributions in the 
fahion 



v( z i,z 2 ) 



/= 



(2/)! 



,-1(1+1) § 



E?lo(2^+l)e- /(/+1)l 



(19) 



It is easy to show that this distribution is normalized to 
unity 



/ 



d 2 Zl d 2 z 2 



p{z\,z 2 ) = 1, 



(20) 



where Z\ and z 2 are given by Eqs. ©, ©, and @2J- 
Note that we must deal with the binomial expression 
(|zi| 2 + |z2| 2 ) 4/ firstly and then integrate over the whole 
complex plane (in two dimensions) in order to verify 
the normalization condition 120JI . The differential ele- 
ment of area in the Z\(z%) plane is d 2 z\ — dxdp x /2H 
(d 2 z 2 — dydp y /2H) 0- Moreover, we have the phase- 
space relationships 
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with 




(21b) 



where a x = a y = ^H/2moj and cr Px = a p — yj hmu/2. 
In Fig. ^ we depict the behavior of the Husimi distribu- 
tion fi(zi,z 2 ) as a function of \z\ at fixed temperature. 
The profile of the Husimi function is similar to Gaussian 
distribution. 



A = 



A={x,p x ,i/,; 



4 



»} 



91n//(zi, z 2 ) 



<9A 



(24) 



where we have introduced a simplified in which the index 
A successively takes the values x,p x ,y, and p y , respec- 
tively. It is easy to prove that the quantity A has the 
following form -see details in appendix Ivl 




where 



r)(zi,z 2 ) 



J=C 



A = r/(zi,z 2 ) 2 



«l 4 



(55-1)1 (2/)! 



(2I)\ 



\ ^ \ 



-/(/+1)0/T 



(25) 



-/(/+i)e/T 



(26) 



Therefore, the corresponding Fisher measure acquires the 
simpler appearance 



j= / —L—l^^^^ rj(zi,z 2 ) 2 , (27) 



FIG. 1: It is depicted the Husimi function fi{z\,Zi) as a 
function of \z\ for several values of the temperature, (T = 
0.16,0,100) for the linear rotator. The behavior of the 
Husimi function looks like Gaussian distribution. The peak 
of the distribution increases as the temperature decreases. 



B. Wehrl entropy and Fisher Information 

The Wehrl entropy is a semiclassical measure of localiza- 
tion 0. So is Fisher's one p|. The Wehrl measure 
is simply a logarithmic Shannon measure built up with 
Husimi distributions. For the present bi-dimcnsional 
model this entropy is of the form 



IV = / d2zid2z % ( Zl ,z 2 )hx^ Zl ,z 2 ), (22) 



where n(z\,z 2 ) is given by Eq. (|19fl . The so-called 
phase-space, shift-invariant Fisher measure \L a is a par- 
ticular instance of the general Fisher-one [!J lld | , that 
can be regarded as an (also semiclassical) counterpart 
of Wehrl entropy [HJ. Extending now to the present 
2D— case the ideas developed in Ref. 13] for the case of 
the ID— harmonic oscillator, we define the (phase-space) 
shift invariant Fisher measure in the fashion 



j= / d2zid2z2 n( Zl ,z 2 )A, (23) 



2 = (77(2:1, z 2 ) 2 ), 
where with the notation 



ir\ f d2zi d2z2 1 \r 

(G) = / p(z)G, 

J ir ir 



(28) 



(29) 



we allude to the semi- classical expectation value of Q . 
In Fig. 12 we plot both Fisher's information and Wehrl's 
entropy as a function of T/0. They behave in different 
manner. If the temperature T — > 0, Fisher's informa- 
tion measure (inverse-delocalization) takes its maximum 
value and Wehrl's its minimum. This behavior is reversed 
for high temperatures, with the degree of delocalization 
becoming larger and larger. 



III. RIGID ROTATOR IN THREE DIMENSIONS 

In the present section we consider a more general prob- 
lem, the 3D— rigid rotator model, whose hamiltonian 
writes Q 



H = 



2U 



f 2 
ZJL 
2L, 



21' 



(30) 



where I x , I yi and I z are the associated moments of in- 
ertia. A complete set of rotor eigenstates is {\IMK)}. 
The following relations apply 
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T/0 

FIG. 2: The trend of the Fisher information measure (de- 
localization) (I) in comparison with the Wehrl entropy (W) 
as a function of the temperature are shown for the linear ro- 
tator. We see that, if the temperature increases, then the 
derealization decreases while Wehrl entropy increases. 



—I, - ■ — 1,0, 1, • • •, I represent the projections on the 
intrinsic rotation axis of the rotor. All states present a 
degeneracy given by (2/+ 1). The parameters I x = I y = 
I xy and I z are the inertia momenta. Several geometrical 
cases are characterized through the I xv /I z ratio. For 
instance, the marginal value I X y/Iz — 1 corresponds to 
the spherical rotor. Limiting cases can be considered; 
this is, Ixy/Iz = 1/2 and I xy /I z — > oo that correspond 
to the extremely oblate and prolate cases, respectively. 

A. Coherent states 

In order to obtain the Husimi distribution for this prob- 
lem we need first of all to construct the associated co- 
herent states, as already discussed by Morales et al. in 
Ref. 0. Introduce first the auxiliary quantity 

X im ,k = + M)\(I - M)\(I + K)\(I - K)\, 

(36) 

and then write 111 



L 2 \IMK) 
L Z \IMK) 
J Z \IMK) 



I(I + 1)\IMK) 

K\IMK) 

M\IMK), 



\z1Z2Z3) = e -2 V LV ; 1 I. 



IMK), (37) 



(3T) where the following supplementary variable were intro- 
duced 



where I = 0, . . . , oo, -I < K < I, and -I < M < 
I. The states \IMK) satisfy orthogonality and closure 
relations 0] 



(I M K \IMK) = 5j> j 5 



m' .m$k',k 



E E E \IMK)(IMK\ = 1. 



1=0 M=-I K=-I 



(32) 



(33) 



If we take L 2 — L 2 +L 2 +L 2 and assume axial symmetry, 
I x — I y , we can recast the hamiltonian as 



i«e«3 I x y 



M 2 = N 2 (1 + M 2 ) 2 (1 + N| 2 ) 2 



(38) 



B. Husimi function, Wehrl entropy, and Fisher 
measure 

Using now Eq. I|37(l we find 

\{IMK\ Zl z,z 3 )\ 2 = ^1^[(2/)!] 2 | 21 | 2 ^ +M '|z 2 | w |^| 2(/+X) 

I .M ,K 

(39) 

and determine that, in this case, the rotational partition 
function reads 



H 



21 



L 



(34) 



where L 2 is the angular momentum operator and L z is 
its projection on the rotation axis z. The concomitant 
spectrum of energy becomes 



£i,k 



21 



Hi + K 



(35) 



where I = 0, 1, 2, • ■ • and it represents the eigenvalue of 
the angular momentum operator L 2 , the numbers m = 



^ = E E E 

1=0 K=-I M=-I 



-Pei,K 



(40) 



Z 3D =£(2/+l) e -^ +1 



1=0 



* t 

K=-I 



.(41) 

It is convenient to remark that if we take one of limiting 
cases, this the extremely prolate, I xy /I z — > oo, the only 
term that survives in the right sum of the right side in 
Eq. (|41|) is that for K — and all terms for K ^ vanish; 
in such particular case Zid is recovered from Z^£>. 
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e-,-1 2 ^ (2/)! 



The pertinent Husimi distribution then becomes 

oo 

E 

where 



ti*u ^,z 3 ) = ^y:^f M 2/ e-V+V9 x g(I) 



z 3D ^ n 

AU 1=0 



(42) 



where a z = ^h/2mu and a Pz = ^Jhmuj /2. In this in- 
stance d 2 z 3 = dzdp z /2h. So, after a bit of algebra we 
arrive to 



i = I dr / i(z 1 ,z 2 ,z 3 )^7 2 (ki| 2 k 2 | 2 + i (i + Ni| 2 ) 2 ) + 



9(1) = £ 7j 
with 



+4M 2 N 2 ^, 



b| 2 = (i + l^l 2 ) 2 !^! 2 



(52) 



i.e.. 



2^2 



h^ = i^(i +\ Z3 n 



(44a) 



(44b) 



I = <7 2 (N 2 M 2 + i (l + N 2 ) 2 ) + 4| U |1z 3 r) (53) 



where 



We can easily verify that (m( z ^i z 2, z 3 ) is normalized in the 
fashion 



-(l + N 2 ) 2 £r=o ^\v\^e-m+^ xg (I) 



£r=o^M 2/ e- /(/+1) 9x 5 (I) 



dT ti^i, z 2 , z 3 ) = 1, (45) 
where dr is the measure of integration given by [| 



dr = dr|4[(l + |z 1 | 2 )(l + |z 3 | 2 )] 4 |z 2 | 4 

|2\l2|_ 2 



" £r=o^M ai e- I(I+1 >*x fl (I) 



(54) 



In the special instance I xy /I z = 1, that corresponds to 
the spherical rotor, we can explicitly obtain 



[(l + |z 1 | 2 )(l + |z 3 | 2 )] 2 |z 2 | 2 + U (46) p(*i>*2>*3)=e 



-\u\ 2 l^i=o i\ K 

E£=o (2/+l) 2 e- / ( / + 1 )f 



(55) 



with 



dr 



d 2 ^i d 2 z 2 d 2 z 3 



7T 7T 7T 



(47) 



Having the Husimi functions the Wehrl entropy is 
straightforwardly computed. 



where, of course, in this case we have three degrees of 
freedom. 

We compute then, the Wherl entropy in the form 



W = / dr^(zi,z 2 ,z 3 ) ln/i(zi,z 2 ,z 3 ), (48) 
and explicitly the Fisher measure as follows 



g 7.5 



1 = J dr>(zi, z 2 , z 3 ) A 3 D 



(49) 




where in this case we define the quantity A 3 d in three 
dimensions as 



A 3D = £ 

^={X-,P^,y,Py,Z,p z } 



(91n/x(zi,z 2 ,z 3 ) 



dA 



, (50) 



with the phase space relationships i)21a.[l . (|21b(l and 



. l2 1 / z 2 pi 
4 \ at at 



(51) 



FIG. 3: A family of Wehrl entropy W as a function of T/Q 
is depicted for the 3D anisotropic rigid rotator for several 
values of the anisotropy (Ixy/Iz = 1/2,1,5,10). The trend 
of the Wehrl entropy is not very similar to the linear case, it 
has a maximum in T/Q = 0, a minimum in T/Q > 2, whose 
exact value depends on the ratio I xy /I z . 

In figure we compare the Wehrl entropy W, as a 
function of T/Q, for several values of I xy /I z ; this is, the 



6 



extremely oblate rotator I xy /I z = 1/2 (e.g., CHCI3 and 
C 6 H 6 ), the prolate rotator I xy /I z — 5,10 (e.g., CH 3 Cl 
and PCI5) and the spherical case I xy /I z = 1. Unfor- 
tunately, from the 3D formulation of coherent states is 
not possible to recover the form of the Wehrl entropy, in 
the same way as occurs with the partition function from 
the Ea. H41|) . which is conveniently lead to the form of 
the Ea. i|18|) in the limiting case of the extremely prolate 
rotator I xy /I z — > 00. It could be due to a note given 
in Ref. 1] , the present version of coherent states formula- 
tion is weak because the measure from the Eq. l|46|l is not 
positive. 

IV. CONCLUDING REMARKS 

We have concentrated our effort on the study of the semi- 
classical behavior of the rigid rotator and have obtained 



in analytical fashion the form of the Husimi distribution 
for two cases, namely, the linear and the axially sym- 
metric rigid rotator. As it is expected, the linear case 
is obtained as a particular case from the formulation in 
three dimensions. Also, we have obtained an analyti- 
cal expression of the shift-invariant Fisher measure built 
up with Husimi distribution, for the rigid rotator model, 
concluding that Fisher measure is better indicator of the 
derealization than Wherl entropy for this model. The 
present study could motivate other specialists to improve 
the present formulation of the coherent states in order to 
recover all quantities in two dimensions from a formula- 
tion in three dimensions. 
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the following result 



51n^(zi,z 2 ) d\z\ 
to ^^to" 



(56) 



where the quantity r](zi,z^) was defined in Eq. 12611 . 
Moreover, from Eqs. (|21a|l and (|21bjl we have 



dx 



(57) 



and we are lead to 



<91n^(zi, z 2 ) _ r\{z\,Z2)x 



dx 



2\z\ai 



(58) 



We arrive to a similar expression differentiating with 
respect to p x , 



91n^(zi,z 2 ) _ ■q(z 1 ,z 2 )p x 



dp x 



2 M< 



(59) 



Analogous expressions are obtained replacing x and p x 
for y and p y . Finally, substituting these results into Eq. 
(|24|l we thus arrive to 



APPENDIX: A BIT OF ALGEBRA 



First of all, we carry out the differentiation the Husimi 
distribution (|19fl with respect to the variable x, obtaining 



A = r](z 1 ,z 2 ) 2 



(60) 



